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The  Counter  Obstacle  Vehicle  (COV)  is  a  heavyweight 
bulldozer-type  vehicle  used  for  earthwork  and  general 
obstacle  removal  in  military  operations.  Microprocessor- 
based  controls  have  been  considered  for  use  on  the  vehicle 
to  automate  blade  control  and  thus  enhance  operator 
performance.  This  study  evaluates  the  feasibility  of  using 
blade  force  measurements  for  controlling  the  depth  of  cut 
in  real-time  operation  of  the  COV.  This  method  of  depth 
control  is  significantly  limited  because  the  variability  of 
natural  soil  deposits  causes  the  blade  force  to  vary  over 
time.  To  account  for  this,  the  study  used  random  process 
(i.e.,  stochastic)  models  of  soil  variation  based  on  empirical 
data  to  simulate  natural  deposits.  Blade  forces  as  the  COV 
plowed  through  the  simulated  soils  were  calculated  using 
soil  mechanics  models.  Kalman  filtering  and  optimal  con¬ 
trol  techniques  were  used  to  adjust  depth  of  cut  based  on 
measured  blade  force.  Study  results  suggest  the  blade  force 
variability  may  not  be  as  great  an  obstacle  to  the  use  of 
force  feedback  for  blade  control  as  previously  thought. 
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ESTIMATING  DEPTH  OF  CUT  FROM  BLADE  FORCES 
ON  THE  ARMY  COUNTER  OBSTACLE  VEHICLE 


1  INTRODUCTION 


Background 

Effective  use  of  the  Counter  Obstacle  Vehicle  (COV)  for  minesweeping  requires 
continuous  plowing  at  a  constant  depth.  However,  because  the  COV  encounters  varia¬ 
tions  in  ground  topography  as  it  moves  along,  the  plow  height  must  be  continuously 
controlled  to  achieve  a  cut  of  constant  depth.  Since  it  is  difficult  for  an  operator  to 
make  such  adjustments  while  closed  inside  the  vehicle,  automatic  control  is  desirable. 

The  first  requirement  for  such  control  is  an  accurate  method  for  monitoring  the 
position  of  the  blade.  For  civilian  vehicles,  front  mounted,  exterior  devices  can  be  used 
to  monitor  blade  position.  As  a  combat  vehicle,  however,  the  COV  operates  in  a 
hazardous  environment  where  devices  mounted  on  the  exterior  of  the  vehicle  are 
vulnerable.  The  combat  environment  also  precludes  the  use  of  visual  inspections  by  the 
COV  operators  or  other  personnel. 

An  attractive  alternative  for  blade  control  is  real-time  force  feedback,  by  which 
the  depth  of  cut  is  estimated  from  the  force  on  the  blade.  However,  the  accuracy  of 
blade  control  using  force  feedback  depends  on  the  accuracy  with  which  plowing  depth  can 
be  estimated  from  blade  force.  Such  estimates  are  complicated  by  the  dependence  of 
blade  force  on  (1)  spatial  variations  of  soil  properties,  (2)  obstacles  such  as  roots  and 
boulders,  and  (3)  dynamic  effects  of  the  blade  passing  through  the  soil. 

Objective 

The  objective  of  this  work  was  to  assess  the  feasibility  of  using  force  feedback  for 
controlling  depth  of  cut  of  the  COV  blade.  Controlling  depth  of  cut  depends  primarily  on 
the  ability  to  accurately  estimate  depth  at  an  arbitrary  time  from  the  time  stream  of 
blade  forces  measured  up  to  and  including  that  time.  Thus,  this  work  addressed  the 
accuracy  of  depth  of  cut  estimates  made  from  blade  forces. 


Approach 

To  evaluate  this  method  of  estimation,  a  model  was  developed  which  is  based  on 
variations  in  the  force  on  the  blade  which  are  caused  by  variations  in  the  properties  and 
topography  of  the  soil.  (The  dynamic  effects  of  the  blade  are  not  modeled.)  To  study  the 
accuracy  of  the  estimate  produced  using  this  model,  several  simulations  were  conducted. 

The  development  of  the  overall  model  begins  with  a  model  of  the  forces  on  the 
blade  (Chapter  2).  These  forces  depend  on  the  characteristics  of  the  soil,  and  especially 
on  the  undrained  shear  strength,  so  that  property  is  modeled  using  stochastic  equations 
(Chapter  3).  Since  the  ultimate  goal  is  to  use  estimates  of  the  force  to  control  the 
position  of  the  blade,  a  model  of  position  and  control  is  developed  using  stochastic 
control  theory  (Chapter  4).  The  final  control  function  will  require  an  estimate  of  the 
blade  position,  which  can  be  obtained  using  the  Kalman  filter  (Chapter  5).  Finally  all 
these  elements,  the  force  and  soil  models,  the  position  equation  from  the  control  model 
(currently  with  zero  external  control),  and  the  Kalman  filter,  are  used  to  simulate  the 
COV  plowing  operation. 
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2  BLADE  LOAD  MODEL 


The  principal  plowing  force  on  the  COV  blade  results  from  shearing  soil  in  front  of 
the  blade.  The  shearing  strength  of  soil  is  usually  modeled  well  by  the  Mohr-Coulomb 
failure  criterion,  expressed  as  1 


Tj.  =  c  +  o^tan(4>) 


[Eq  11 


where  ij 


shear  stress  on  the  failure  surface 


Oj  =  normal  stress  on  the  failure  surface 
c  =  undrained  soil  strength  (cohesion)  and 
=  soil  friction  angle. 

The  soil  parameters  c  and  <t>  usually  vary  in  space. 

In  pushing  the  blade  into  the  soil,  horizontal  stresses  in  the  soil  mass  increase  until 
a  failure  surface  develops.  On  this  surface  the  stresses  x  and  o  satisfy  the  Mohr- 
Coulomb  criterion  (Figure  1).*  The  mechanics  of  blade  action  and  soil  failure  in  this 
pushing  case  are  similar  to  passive  failure  of  retaining  structures.  Considering  two- 
dimensional  equilibrium  of  the  failed  mass  of  soil,  the  force  exerted  on  the  blade  can  be 
expressed  as2 


F  =  0.5 Yd  N  +  2cdN  +  qdN  [Eq  2] 

(|)  C  CJ  ~ 

where  F  =  force  on  the  blade 
y  =  unit  weight  of  soil 
d  =  depth  of  cut 
c  =  cohesion  intercept 
q  =  surcharge  on  the  soil  surface. 

,  Nc,  and  N  are  constants  which  depend  on  the  soil  friction  angle  and  the  blade's 
roughness.  This  force  model  is  the  sum  of  three  terms.  The  first  term  is  the  contribu¬ 
tion  of  lifting  the  failed  mass  up  the  failure  surface,  the  second  term  is  the  undrained 
shearing  of  the  soil,  and  the  third  term  is  the  effect  of  surcharge. 

Equation  2  is  simple  in  form,  but  the  variability  of  the  soil  parameters  y,  c, 
N^,  Nc,  and  make  the  estimation  of  depth  from  force  measurements  impractical. 
Simplifications  are  necessary  to  apply  this  equation  for  the  COV  application. 

Simplifying  Assumptions 
Smooth  Blade 

A  smooth  blade  means  that  the  soil  exerts  no  lateral  traction  on  the  blade.  Soil 
will  not  stick  to  the  blade  as  it  is  plowed.  Assuming  low  friction  between  the  blade  and 
soil,  the  parameters  N  ,  N  ,  and  N  simplify  to  functions  of  the  soil  friction  angle. 


‘T.  W.  Lambe  and  V.  R.  Whitman,  Soil  Mechanics  (John  Wiley  <5c  Sons,  New  York,  1968). 
*  All  tables  and  figures  can  be  found  at  the  end  of  this  report  (p  27). 

2Lambe  and  Whitman. 
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N  =  tan2(45°  +  <t> / 2 )  [Eq  3] 

4> 

N  =  N  =  tan(45°  +  / 2 )  [Eq  4] 

c  q 

and  the  failure  surface  is  a  plane  inclined  at  (45°  +  4>/2  )  from  the  vertical  blade  surface. 
Undrained  Condition 

"Undrained  condition"  means  that  the  soil  is  plowed  quickly  enough  that  water  is 
not  squeezed  out  of  the  soil  during  shearing.  Thus,  the  water  content  of  the  soil  is 
constant  during  plowing.  If  it  is  assumed  that  the  soil  is  sheared  at  constant  water 
content  and  that  the  blade  is  smooth,  then  Equation  1  reduces  to  t  =  c,  and  the  blade 
load  (Equation  2)  becomes 

F  =  0.5yd2  +  2cd  +  qd  [Eq  5] 

where  the  cohesion  intercept,  c,  is  the  undra  ned  shear  strength.  The  assumption  of 
undrained  shearing  of  soil  is  reasonable  for  the  COV  operating  at  design  plowing  speeds 
or  in  low  permeability  soils  such  as  silts  and  clays. 

To  examine  the  importance  of  each  term  in  the  force  model  expressed  in  Equation 
2,  a  typical  case  was  chosen  with  a  blade  depth  of  2  ft*  and  loose  soil  in  front  of  the 
blade.  The  geometry  of  this  case  is  shown  in  Figure  2.  Choosing  typical  c/y  ratios  in 
the  range  2  to  6,  the  contribution  of  each  term  as  a  percentage  of  total  force  is  shown  in 
Table  1. 

At  average  strength  values  (c/y  =  4  to  5)  the  dominant  term  in  the  force 
function  is  the  undrained  shearing  (2cd).  As  the  depth  increases,  the  contribution  of  the 
first  term  (0.5yd2)  increases,  but  the  limit  on  blade  height  places  a  physical  limit  on 
the  magnitude  of  this  term.  For  blade  heights  up  to  4  ft,  shearing  is  the  dominant 
component  of  force.  This  observation  can  be  interpreted  to  mean  that  for  optimal 
control,  variability  in  unit  weight  y  can  be  ignored  compared  to  variability  in  shear 
strength:  at  typical  operating  depths  (<  2  ft),  random  variations  in  unit  weight  result  in 
insignificant  changes  of  force.  For  example,  at  c/y  =  4,  a  10  percent  change  in  unit 
weight  causes  a  1  percent  change  in  force. 

Comparisons  With  Field  Measurements 

To  confirm  a  force  model  similar  to  Equation  5  (with  q  =  0),  O'Callaghan  and 
McCullen  performed  full-scale  bulldozer  tests.3  A  dynamometer  was  mounted  on  the 
bulldozer  blade  and  the  force  exerted  on  the  blade  was  measured  for  various  depths. 
Measured  and  calculated  values  of  horizontal  blade  force  are  compared  in  Figure  3.  The 
theoretical  relationship  is  approximately  linear,  and  thus  satisfactory  for  present  pur- 

2 

poses,  because  of  the  insignificant  contribution  of  the  lift  term  (0.5yd  ). 


*A  metric  conversion  table  is  provided  at  the  end  of  this  report  (p  25). 

3J.  R.  O'Callaghan  and  P.  J.  McCullen,  "Soil  Mechanics  in  Relation  to  Earth-Moving 
Machinery,"  Institution  of  Mechanical  Engineers— Proceedings  (Symposium  on  Earth- 
Moving  Machinery),  Vol  179,  part  3F  (1964-65),  pp  23-29. 
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3  MODELING  SOIL  VARIABILITY 


Covariance  Function  for  Undrained  Shear  Strength 

In  the  force  model  outlined  in  Chapter  2,  the  effects  of  soil  are  represented  by  y, 
the  unit  weight  and  c,  the  undrained  shear  strength.  The  latter  quantity,  undrained 
shear  strength  of  the  soil,  varies  horizontally  and  vertically  in  any  deposit,  even  those 
said  to  be  homogeneous.  This  variability  is  usually  modeled  as  a  second-order  stochastic 
field  with  statistical  parameters  estimated  from  field  data.4  A  second-order  stochastic 
field  is  similar  to  models  used  in  other  branches  of  science  to  represent  fluctuating  time- 
series  of  data  (e.g.,  economics),  only  a  stochastic  field  is  defined  in  space  rather  than 
time  and  can  be  defined  over  multiple  dimensions.  Stochastic  field  models  are  used  to 
characterize  the  spatial  variability  of  soils  because  no  amount  of  field  testing  can  ever 
completely  characterize  the  point-to-point  variation  of  natural  soils. 

A  second-order,  stationary,  random  process  is  completely  specified  by  a  mean  (or 
mean  trend)  and  a  covariance  function.  The  function  describes  the  degree  of  covariance 
between  soil  properties  as  a  function  of  separation  in  space.  For  a  stationary  process, 
the  covariance  between  any  two  points  depends  only  on  the  (vector)  separation  between 
points  and  not  on  their  location  within  the  domain.  Thus,  the  covariance  function  is  the 
same  everywhere.  Experience  with  undrained  shear  properties  of  natural  soils  has 
demonstrated  that  the  covariance  function  can  be  adequately  modeled  as  an  exponential 
function  of  the  form 

B(r)  =  o^e  r/ro  [Eq  6] 


where  B(r) 


r 


covariance  function  evaluated  at  separation  r 
variance  of  the  fluctuations  around  the  mean 
a  scale  parameter  called  the  covariance  distance 
separation  between  any  two  points. 


Empirically,  the  covariance  distance  rQ  is  usually  anisotropic,  shorter  by  a  factor  of 
ab^ut  10  in  the  vertical  compared  to  horizontal  direction. 

De  Groot  used  a  maximum  likelihood  method  to  estimate  the  covariance  parame- 
2 

ters,  a  and  r  ,  for  several  natural  soil  deposits.  His  results  suggest  that  horizontal 
o 

covariance  distances  in  different  soil  deposits  can  vary  from  50  to  200  ft  and  that 
coefficients  of  variation  (standard  deviation  divided  by  mean)  range  from  15  to  35 
percent.  These  results  are  consistent  with  the  findings  of  others.5 


“D.  J.  DeGroot,  "Maximum  Likelihood  Estimation  of  Spatially  Correlated  Soil  Proper¬ 
ties,"  M.S.  Thesis,  M.I.T.,  May  1985. 

5G.  B.  Baecher,  et  al.,  "Geotechnical  Reliability  of  Offshore  Gravity  Platforms,"  M.I.T. 
Sea  Grant,  Report  MITSG  80-20,  Cambridge,  MA,  1980;  Spikula,  D.,  "Statistical 
Estimation  of  Soil  Properties:  An  Application,"  M.S.  Thesis,  M.I.T.,  May  1983;  W.  Tang, 
"Probabilistic  Evaluation  of  Penetration  Resistance,"  American  Society  of  Civil 
Engineers  Journal,  Geotechnical  Engineering  Division,  Vol  105,  GT10  (1979)  pp  1173- 
1191. 
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Covariance  functions  and  coefficients  of  variation  reported  in  the  literature  usually 
refer  to  point  properties  measured  in  individual  tests  or  at  specific  depths  in  a  boring. 
As  such,  they  summarize  the  stochastic  variability  among  small  elements  of  soil.  A 
blade,  on  the  other  hand,  moves  soil  within  a  much  larger  volume.  Variations  in  blade 
load  are  related  to  variations  among  the  average  properties  of  such  larger  soil  volumes. 
Since  variations  from  point  to  point  within  these  larger  volumes  are  averaged,  the 
coefficient  of  variation  of  volume-averaged  properties  is  reduced  compared  to  those  of 
point  properties.  At  the  same  time,  the  covariance  function  is  spread  out,  giving  larger 
covariance  distances. 

The  reduction  of  variability  with  volume  and  the  spreading  of  the  covariance 
function  is  counteracted  by  discrete  changes  in  soil  type  that  may  occur  along  the  path 
of  plowing.  Soil  properties  across  such  discrete  changes  have  little  correlation  with  one 
another,  and  thus  both  increase  the  volume-averaged  coefficient  of  variation  and 
decrease  the  covariance  distance. 

In  the  simulation  studies  it  was  decided  to  strike  a  balance  between  these  opposing 
trends  by  adopting  a  coefficient  of  variation  of  10  percent  and  a  covariance  distance,  r  , 
of  50  ft.  Choosing  a  typical  average  undrained  shear  strength  of  400  psf,  the  stochastic 
behavior  of  the  undrained  shear  strength  can  be  completely  specified  by  its  first  and 
second  moments: 


where  E 
r 


E[c]  =  400  psf 

B  (r)  =  1600e“r/50 
cc 

the  expectation  (i.e.,  mean  or  arithmetical  average) 
the  separation  in  feet. 


[Eq  7] 
[Eq  8] 


Stochastic  Differential  and  Difference  Equations  for  Undrained  Shear  Strength 

Optimal  stochastic  control  (discussed  in  Chapter  4)  requires  that  random  processes 
must  be  specified  in  stochastic  differential  or  difference  form.  One  of  the  many 
stochastic  processes  of  this  form  is  the  first-order  Gauss-Markov  process,  which  is  a 
normal  or  Gaussian  process  with  an  exponential  covariance  function.  The  covariance 
function  of  the  undrained  shear  strength  of  soil  (Equation  6)  is  also  exponential,  so  it  is 
assumed  that  the  undrained  shear  strength,  like  the  Gauss-Markov,  is  a  normal  process. 
Therefore,  its  behavior  can  be  described  by  a  stochastic  differential  equation  for  a 
normal  process:6 

dc/dt  =  -(l/ro)c  +  (L/ro)E[c)  +  e  [Eq  9] 


where  c  =  undrained  shear  strength  of  the  soil 

rQ  =  covariance  distance  of  the  undrained  shear  strength  of  the  soil. 


6A.  Gelb,  Applied  Optimal  Estimation  (M.I.T.  Press,  Cambridge,  MA,  1974). 


In  Equation  9,  e  is  Gaussian  white  noise  with 

E[e]  =  0 

E[e2]  =  (2/r  )o2 
o  c 

where  E  =  the  expectation  operator 

rQ  =  the  covariance  distance  of  the  undrained  shear  strength  of  the  soil 

a  „  =  the  variance. 


# 


In  discrete  time  Equation  9  can  be  written  as  a  stochastic  difference  equation, 

c(t+l)  =  yc(t)  +  v(t)  [Eq 

where  v(t)  is  an  independent  Gaussian  sequence,  and  from  Gelb  and  Equations  7  and  8 


T  =  e  1/ro  =  0.98 

E[v2]  =  a2  (1  -  Y2)  =  62.7. 
c 


4  BASICS  OF  STOCHASTIC  CONTROL  THEORY 


Stochastic  control  theory  deals  with  dynamic  systems  whose  states  are  described  by 
stochastic  differential  or  difference  equations.  For  example,  simple  stochastic  control 
involves  a  linear  system  disturbed  by  Gaussian  white  noise.  Linearity  refers  to  measure¬ 
ments,  y(t),  of  the  system  state,  x(t),  which  depend  linearly  on  the  values  of  x(t). 
White  noise  refers  to  variations  in  the  state  variables,  x(t),  which  are  independent  and 
identically  distributed  for  different  values  of  t.  The  objective  of  this  theory  below  is  to 
determine  a  function  [u(t)  below]  that  will  "control"  the  state  of  the  system  [described 
by  x(t)]  so  as  to  minimize  the  deviations  of  that  state  from  some  ideal  value  or  set  of 
values.  There  are  two  cases  to  be  considered:  the  ideal  case,  where  measurements  of 
the  state  of  the  system,  x(t),  are  exact,  thus  the  only  noise  is  the  Gaussian  white 
noise,  w(t);  and  the  "noisy"  case,  where  observations  of  the  system,  j(t)  ,  contain 
noise  v(t),  in  addition  to  white  noise.  The  derivations  below  address  both  cases. 


System  with  White  Noise  Only 

In  the  first  case  it  is  presumed  that  the  initial  conditions  are  random,  but  that  the 
state  of  the  system  can  be  measured  with  complete  certainty  at  any  time.  The  system 
can  be  represented  by  the  equation 

x(t+l)  =  yx(t)  +  Tu(t)  +  w(t)  [Eq  15] 

where  x  =  n-dimensional  state  vector 

u  =  m-dimensional  control  vector 
w  =  n-dimensional  vector  of  random  perturbations 
y  =  n*n  one-step  transition  matrix  from  x(t)  to  x(t+l) 
r  =  n*m  one-step  transfer  matrix  from  u(t)  to  x(t+l). 

By  the  Gaussian  white  noise  assumption,  the  random  vectors  w(0),  w(  1 ),.... ,  w(t)  are 
independent,  with  mean  E[w(t)]  and  covariance  E[w(t)  wT(t)]  having  values, 

E[w(t ) ]  =  0,  E[w(t )  wT(t)]  =  Q  [Eq  16] 

where  T  is  the  transpose.  The  initial  condition  x(0)  is  taken  to  be  unknown,  but  with 
specified  mean  E[x(0))  and  covariance  matrix  E[x(0)xT(0)]  having  the  values, 

E[x(0) ]  =0,  E[x(0)  xT(0)]  =  Pq  [Eq  17] 


The  control  procedure  attempts  to  minimize  deviations  of  the  process  x(t)  from 
some  ideal  value  or  trajectory  of  values.  Without  loss  of  generality,  these  ideal  values 
are  set  to  zero.  For  practical  reasons,  squared  deviations  of  x(t)  from  the  ideal  values 
x(t)  =  0  are  frequently  chosen  to  be  minimized.  This  leads  to  a  performance  criterion, 
J,  of  the  quadratic  form, 

J  =  e{£ [x^( t )Ax(t )  +  u^(t)Bu(t)]}  [Eq  18] 

where  A  and  B  are  positive  definite  matrices  which  express  the  cost  of  deviations  of  the 
components  of  x(t),  and  the  cost  of  controls,  respectively.  Mathematically,  the 


matrices  A  and  B  are  arbitrary.  Quantitative  values  for  A  and  B  depend  on  the  applica¬ 
tion  in  question. 

Because  this  "ideal"  situation  is  identical  to  a  deterministic  control  problem,  the 
controller  u(t)  which  minimizes  J  (Equation  18)  is7 


u( t )  =  -Lx(t) 

where  L,  the  feedback  matrix,  depends  on  ¥,  r,  A,  and  B. 


[Eq  19] 


System  with  Additive  Noise 

For  most  optimal  control  problems,  including  force  feedback  for  the  COV,  the  state 
variables,  x(t),  are  observed  indirectly  and  with  measurement  errors.  Let  ^(t)  be  the 
vector  of  observations  of  the  states  x(t);  then 

^(t)  =  g(x(t))  +  v(t)  [Eq  20] 

where  g(.)  =  a  deterministic  function 

v(t)  =  a  vector  of  independent,  identically  distributed  (iid)  measurement 
errors. 

In  the  blade  control  problem,  the  function  g(x(t))  is  the  blade  load  model  of  Equation  5, 
or  a  linearized  version  of  it.  The  state  vector  comprises  blade  depth  and  undrained  shear 
strength.  In  most  cases  j>(t)  is  a  linear  function  (discussed  in  the  next  section)  such  that 
Equations  15  and  20  lead  to  the  system  of  equations, 


x(t+l)  =  fx(t)  +  ru(t)  +  w(  t ) 
^(t)  =  Hx( t )  +  v(t) 


[Eq  21] 


where  H  is  the  deterministic  matrix  of  coefficient  terms  expressing  the  linear  function 
g(.),  and  w(t)  and  v(t)  are  independent  Gaussian  white  noise  sequences  with  moments 


E( w( t ) ]  =  0,  E[w(t)  w  ( c ) ]  =  Q 

E[v(t)]  =  0,  E[v(t)  vT(t)]  =  R 


[Eq  22] 


and  x(0 )  is  a  Gaussian  vector,  independent  of  v(t)  and  w(t),  with  properties  defined 
in  Equation  17. 


7 A.  E.  Bryson  and  Y.  C.  Ho,  Applied  Optimal  Control  (Blaisdell  Publishing  Co.,  Waltham, 
MA,  1969). 
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The  optimal  controller  for  the  system  of  Equation  21,  given  the  performance 
criterion  of  Equation  18,  is8 

u(t)  -  -  Lift)  [Eq  23] 

where  L  =  the  feedback  matrix 

x(t)  =  the  mean  estimate  of  x(t)  conditioned  on  the  measurements 
y(0) ,  y( 1 ) , . . . ,y( t ) . 

The  estimate  x(t)  is  used  as  if  it  were  exact,  similar  to  deterministic  control. 

The  solution  of  Equation  23  results  from  the  so-called  "separation  theorem."9  This 
theorem  holds  that  stochastic  control  problems  with  a  quadratic  performance  criterion 
(e.g.,  Equation  18)  can  be  separated  into  two  independent  parts:  one  is  an  optimal  filter 
by  which  x(t)  is  estimated  from  the  observations  ^(t);  the  other  is  linear  feedback 
from  the  estimated  state,  x(t),  obtained  in  the  same  way  as  in  the  case  of  no  observa¬ 
tion  errors. 

The  process  by  which  x(t)  is  estimated  is  based  on  Kalman  filtering.  The  filter  is 
recursive  and  depends  only  on  the  system  dynamics  and  noise  characterization.  The 
estimate  x(t)  is  independent  of  the  quadratic  performance  criterion  J.  Thus,  controlled 
behavior  for  a  dynamic  system  depends  to  a  large  extent  on  the  performance  of  the 
Kalman  filter.  As  a  result,  a  major  portion  of  this  feasibility  analysis  has  concentrated 
on  the  adequacy  of  Kalman  filter  estimation  of  depth  from  the  time  stream  of  the  blade 
load. 


Nonlinear  Observations 

In  many  physical  systems,  including  the  blade  load  model  of  Equation  5,  the 
observation  vector  y(t)  is  a  nonlinear  function  of  the  state  variables  x(t).  To  obtain 
practical  solutions,  the  nonlinear  function  £(.)  can  be  made  linear  by  expanding  it  in  a 
Taylor  series  around  a  known  vector,  for  example,  the  expected  value  of  x(t).10 
Filtering  and  control  operations  are  performed  on  this  linearized  system.  This  approach 
is  used  in  Chapter  6,  where  the  blade  load  model  is  linearized  around  the  desired  depth  of 
cut  and  the  mean  undrained  shear  strength. 


8K.  J.  Astrom,  Introduction  to  Stochastic  Control  Theory  (Academic  Press,  New  York, 
NY,  1970). 

9K.  J.  Astrom. 

10  A.  Gelb. 


5  KALMAN  FILTER  -  THEORY  AND  RESULTS 


This  section  summarizes  rudimentary  Kalman  filter  theory  as  applied  to  the 
estimation  of  blade  depth  from  force  measurements.  The  Kalman  filter  is  used  to 
estimate  the  conditional  mean  of  the  state  variables,  x(t),  from  which  the  controller  is 
calculated  as  in  Equation  23.  Clearly,  the  accuracy  of  the  control  signal  u(t)  and 
the  behavior  of  the  controlled  system  x(t)  depend  on  the  accuracy  with  which  x  ( t  )  is 
estimated. 

Force  on  the  blade  is  a  function  of  two  random  variables,  the  undrained  shear 
strength  of  the  soil  being  plowed  and  the  depth  of  cut.  The  greater  the  soil  strength  and 
the  deeper  the  cut,  the  greater  the  force.  The  force  feedback  system  simultaneously 
estimates  two  unknowns  and  makes  adjustments  to  (i.e.,  controls)  the  blade  accordingly. 
The  extent  to  which  the  depth  of  cut  at  time  t+1  can  be  precisely  controlled  depends  on 
the  capability  of  the  filter  to  adequately  estimate  the  existing  depth  of  cut  at  time  t. 
This  estimate  is  hindered  by  random  variations  of  the  undrained  shear  strength  and 
inaccuracies  in  the  model  g( . ) . 


Kalman  Filter 

Kalman  filtering  arises  in  many  engineering  problems  where  an  estimation 
algorithm  is  required  to  sequentially  process  observations  of  a  dynamic  system.  Since 
their  development  in  the  early  1960s,1  1  the  properties  of  the  Kalman  filters  have  been 
well  studied. 

From  Equation  15  the  process  equation  for  an  uncontrolled  system  is 

x(t  +  l)  =  't'x(t)  +  w(t)  [Eq  24] 

where  x(t)  =  the  state  vector  at  time  t 
f  =  a  transition  matrix 

w(t)  =  a  vector  of  iid  perturbations  (i.e.,  white  noise). 

The  state  x(t)  is  observed  by  noisy  measurements  (the  observation  includes  additive 


noise): 

^( t )  =  Hx( t )  +  v( t ) 

[Eq  25] 

where 

j[(t)  =  the  observation  vector 

H  =  a  deterministic  matrix 

v(t)  =  a  vector  of  iid  observation  errors. 

11 R.  E.  Kalman,  "A  New  Approach  to  Linear  Filtering  and  Prediction  Problems," 
American  Society  of  Mechanical  Engineers  (ASME)  Journal  of  Basic  Engineering ,  Vol 
82  (1960),  pp  34-45;  R.  E.  Kalman  and  R.  S.  Bucy,  "New  Results  in  Linear  Filtering  and 
Prediction  Theory,"  ASM E  Journal  of  Basic  Engineering,  Vol  83  (1961),  pp  95-107. 
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The  initial  state  x(0),  perturbations  w(t),  and  measurement  errors  v(t)  are  assumed 
to  be  uncorrelated,  with  moments 


E[x(0) ]  =  0, 
E[w(t ) ]  =  0, 
E[w( t ) ]  =  0, 


E[x(0)x(0)  ]  =  P 


E(w(t)w(t)  J  =  Q 


[Eq  26] 


E[v(t)v(tr]  =  R 


The  Kalman  filter  can  be  obtained  as  the  linear  Bayes'  least  squares  estimate  of  the 

state  vector  x(t),  given  the  observations  ^(0),  ^(1), . ,£(t).  12  To  indicate 

the  time  at  which  an  estimate  of  x(t)  is  made— because  different  amounts  of  informa¬ 
tion  are  available  at  different  times— let  the  estimate  conditional  on  observations  up  to 
time  t  be  denoted  x(t/t),  rather  than  x(t)  as  used  in  the  previous  section.  In  this 
manner,  an  estimate  conditioned  on  observations  up  to  time  (t-1)  would  be  denoted 
x(t/t-l).  The  equations  of  the  Kalman  filter  can  be  written  as 


x(t/t)  =  x(t/t-l)  +  K(t)[y(t)  -  Hx( t/t-l ) ] 


[Eq  27] 


where  x(t/t-l)  =  estimate  of  x(t),  given  y(t-l) 

=  ¥x( t-l/t-1 ) 

=  Kalman  gain  matrix 
=  P(t/t-l)HT{H[P(t/t-l)]HT  +  R}-1 


[Eq  28] 


K(t) 


[Eq  29] 


where  H 


=  matrix  relating  the  measurements  y(t)  to  the 
measured  state  variables  x(t), 


P(t/t-l) 


=  covariance  matrix  of  the  estimate  x( t/t-l)  in 
Equation  28 


=  [p(t-i/t-i)],r  +  q 


[Eq  30] 


P(t/t) 


covariance  matrix  of  the  estimate  x(t/t) 
in  Equation  27 


[Eq  31] 


(I  -  K(t)H)[P(t/t-l)] 


where  I  is  the  identity  matrix. 


The  Kalman  filter  is  recursive  in  that  the  estimate  of  the  state  variable  at  time  t, 
x(t/t),  is  completely  defined  from  the  observation  vector  ^(0  and  the  previous  step 
estimate  x( t-1  /t-1).  At  time  zero,  the  covariance  matrix  of  x(0)  is  given  by  Equation 
26.  This  is  used  to  calculate  the  covariance  matrix  of  the  first  time  step  estimate  before 
a  measurement  is  made,  P{l/0),  using  Equation  30.  The  covariance  matrix  after  the  first 
measurement  is  made,  P(l/1),  is  calculated  using  Equation  31.  Subsequent  covariance 
matrices  up  to  P(t/t)  are  calculated  by  recursive  application  of  Equations  30  and  31. 


1 2 A.  Gelb. 
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This  recursive  property  reduces  computation  and  computer  memory  requirements. 
Also,  the  estimate  variances  and  covariances  of  Equations  29  through  31  are  independent 
of  the  estimates  in  Equations  27  and  28.  This  means  that  covariance  matrices  and  filter 
gains  can  be  calculated  in  advance  and  stored  in  the  computer  system.  Thus,  the 
recursive  calculations  do  not  have  to  be  made  while  the  system  is  in  operation  (i.e.,  while 
the  COV  is  moving),  thus  real-time  application  of  the  Kalman  filter  is  enhanced. 

Kalman  Filter  Applied  to  a  Linear,  Scalar  System 

To  illustrate  the  performance  of  Kalman  filters  for  simple  system  dynamics,  a 
series  of  simulations  were  performed.  The  system  model  was  taken  as  linear  and  scalar: 


x(t+t)  =  yx(t)  +  w(t) 
^(t)  =  x(t)  +  v(t) 

where  the  notation  follows  that  of  Equations  15  and  20. 


Let  the  covariance  function  of  x  be 

B  (r)  =  a2  e  r/ro 

XX  X 

where  Bxx(r)  =  the  covariance  of  x(i),  x(i+r) 

ox  =  the  variance  of  x  across  its  domain 

r  =  the  separation  distance 

rQ  =  the  covariance  distance  (a  soil  parameter). 

For  the  covariance  function  of  Equation  33,  the  system  components  become 


u,  -  -r/r 
V  =  e  o 


E[ w2 ( t ) ]  =  Q  =  o2(l  -  y2) 


E [ x ( 0 ) )  =  0,  E[x2(0) ]  =  P  =  o2 

O  X 


E( v  (t)]  =  R  =  constant. 


[Eq  32] 


[Eq  33] 


[Eq  34] 


[Eq  35] 


[Eq  36] 


[Eq  37] 


Table  2  is  a  summary  of  seven  simulated  cases.  Case  3  is  designated  the  base  case 
because  its  covariance  distance  is  that  used  for  the  depth-of-cut  model  in  later 
sections.  The  other  cases  show  sensitivity  of  filter  estimates  to  covariance  distance  rQ, 
observation  variance  R,  and  variance  of  the  state  o2.  The  last  column  (designated  P  ) 
is  the  estimation  variance,  P(t/t),  for  large  time.  P^  is  a  measure  of  the  limiting 
precision  of  the  estimate.  Figures  4  through  10  show  estimated  vs.  true  values  of  x  and 
time  plots  of  those  values. 
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Figures  4  through  6  show  sensitivity  with  respect  to  the  covariance  distance,  and  as 
a  result  the  perturbation  variance  Q;  i.e.,  as  the  covarianee  distance  (rQ)  increases,  the 
estimates  increase  in  accuracy  (P  decreases).  This  is  caused  by  the  higher  signal-to- 
noise  ratio  (  Pq/Q  )  as  Q  decreases.  In  the  limiting  case  rQ+  ®  and  the  process  becomes 
a  random  constant  with  P  *0  (after  a  large  number  of  observations  the  constant  is 

co 

exactly  estimated). 


Cases  6  and  7  demonstrate  the  analogous 
(P  =  ah  decreases  the  accuracy  of  the 
The  implications  of  these  mathematical  properties  for  blade  control  are 


As  expected,  cases  3,  4,  and  5  show  that  an  increase  in  observation  variance,  R, 
increases  the  uncertainty  in  the  estimate  of  x. 
result  that  an  increase  in  the  fluctuation  of  x 
estimate. 

consistent  with  the  intuitive  notions  that  (1)  the  less  precisely  blade  force  are  measured, 
the  less  precisely  depth  of  cut  can  be  estimated,  and  (2)  the  more  variable  soil  properties 
are,  the  less  precisely  depth  of  cut  can  be  estimated.  The  behavior  of  the  asymptotic 
estimator  variance  P  indicates  the  performance  properties  of  the  estimator  when  the 

CD 

COV  has  plowed  a  long  way  (i.e.,  when  initial  effects  have  been  attenuated). 


6  ESTIMATING  DEPTH  OF  CUT  FROM  BLADE  LOAD 


This  chapter  presents  results  of  simulation  studies  to  assess  the  feasibility  of 
estimating  the  depth  of  cut  from  measured  blade  force.  The  problem  is  modeled  as  a 
stochastic  control  (estimation)  problem,  where  the  state  variables  are  depth  of  cut  and 
undrained  shear  strength  of  the  soil,  and  force  is  a  nonlinear  response  to  the  state 
variables. 


System  Model 

The  dynamics  of  depth  of  cut  depend  on  mechanical  characteristics  of  the  blade, 
the  dynamics  of  control,  and  the  magnitude  of  noise  in  the  system.  The  difference 
equation  governing  depth  of  cut  is  taken  as 

d( t  +  1 )  =  f d<l( t )  +  fu( t )  +  «d(t)  [Eq  38] 

where  d(t)  =  the  (scalar)  depth  of  cut  at  time  t 
u(t)  =  the  (scalar)  control  input 
wrf(t)  =  (scalar)  random  noise. 

Note  that  Equation  38  implies  that  the  depth  of  cut,  at  zero  control  [wd(t)  =  Oj, 
follows  a  one-step  Gauss- Markov  process. 

In  reality  the  dynamics  of  the  blade  may  differ  from  the  system  of  Equation  38. 
However,  as  long  as  the  state  equation  can  be  presented  in  a  linear  stochastic  difference 
equation,  the  details  of  the  estimation  process  are  only  slightly  changed. 

For  the  present  simulation  studies,  the  following  parameters  were  taken  as 
representative. 

E[d(  t )  1  =  2  ft,  F.[(d(t)  -  2 ) 2  ]  =  P  =  0.50  ft2 

o 

fd  =  e"l/ro  =  e“1/10  =  0.905  [Eq  39] 

E [ w  , ( t ) ]  =  0,  E[w2  ( t ) ]  =  (1  -  Y2)  =  P  =  0.091 
a  a  do 

Note  that  the  chosen  parameters  indicate  that  the  depth-of-cut  process  has  a  mean  of 
2  ft,  covariance  distance  of  10  ft,  and  standard  deviation  of  0.5  ft.  The  mean  (2  ft)  is 
approximately  the  design  specification  (target  value)  for  mine  plowing,  the  standard 
deviation  (0.5  ft)  is  a  conservative  guess  at  the  uncertainty  in  initial  depth  of  cut,  and 
the  covariance  distance  (10  ft)  is  a  conservative  estimate  presuming  on  erratic  topo¬ 
graphy. 
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In  the  simulation  study  it  was  assumed  that  the  control  signal,  u(t)  is  zero, 
indicating  a  fixed  setting  for  the  blade.  For  this  case,  random  variations  of  blade  depth 
are  due  to  noise  and  changes  in  topography.  The  addition  of  control  terms  u(t)  is  to  be 
part  of  a  second  phase  of  the  work  reported  here. 


Soil  Model 

Several  assumptions  were  made  to  model  the  soil.  A  balance  was  struck  between 
the  reduction  of  variability  with  volume  and  discrete  changes  in  soil  type  (see  Chapter  3) 
by  assuming  a  coefficient  of  variation  of  10  percent,  a  covariance  distance,  rQ,  of  50  ft, 
and  an  average  undrained  shear  strength  of  400  psf. 


Linearization  of  Force  Model 

The  force  function  which  was  developed  in  Chapter  2, 

F  =  0.5yd2  +  2cd  +  qd  [Eq  40] 

in  a  nonlinear  function  of  the  state  variables  d  =  depth  of  cut  and  c  =  undrained  soil 
strength  (cohesion).  This  is  the  function  used  to  relate  observations  ^(t)  to  the  state 
variables  x(t).  To  make  the  force  function  linear,  it  was  expanded  in  a  Taylor  series 
about  the  mean  values  of  d  and  c  and  truncated  to  the  first  two  terms.  Denoting  the 
mean  values  of  depth  and  undrained  shear  strength  by  md  and  , 

F  =  F(m  ,m  )  +  (d  -  m  )-f|  +  (c  -  m  [Eq  41] 

d  c  d  od  C  oc 

where  the  partial  derivatives  are  evaluated  at  m  and  m, .  Choosing  y  =120  psf, 

c  d 

q  =  0  (no  surcharge),  and  the  mean  values  m,  =2  ft  and  m  =  400  psf,  the  specific 

d  c 

blade  force  (lb/unit  width)  is, 

F  *  -1840  +  1040d  +  4c  [Eq  42] 

To  assess  the  magnitude  of  error  in  the  linearized  form,  several  comparisons  were 
made.  Figure  11  shows  the  comparison  for  c  =  400  psf,  the  average  soil  strength.  The 
approximation  is  good  for  depths  in  the  range  1  to  3  ft.  Figure  12  is  a  comparison  of 
estimated  vs.  true  depth  in  which  depth  is  back-calculated  from  the  linearized  force, 
d  =  (F  +  240)/1040.  Knowing  the  true  depth,  a  force  is  calculated  from  Equation  5,  then 
the  linear  approximation  is  used  to  back-calculate  an  approximate  depth;  then  real  and 
approximated  depths  are  compared.  The  back-calculated  depth  in  Figure  12  is  close  to 
the  true  depth  in  the  range  0.8  to  3  ft. 

These  comparisons  indicate  that  the  error  in  the  linearized  force  function  is 
small.  However,  more  robust  linearization  techniques  are  available  if  needed. 1 3 


Simulation  Case  Studies 


Based  on  the  above  discussion,  the  following  system  and  observation  equations  were 
used  to  characterize  the  estimation  problems. 


d(t  +  l )  =  ¥,d(t)  +  w,(t) 
a  a 

[Eq  43] 

c(t+l )  =  ¥  c(t)  +  w  (t) 
c  c 

[Eq  44] 

•1840  +  1040d(t )  +  4c(t)  +  v_(t) 

F 

[Eq  45] 

Parameter  values  are  taken  from  Equations  39  and  42.  Equation  45  is  the  observation  of 
the  state  parameters  through  the  force  model,  where  a  random  noise  term  Vp(c)  ex¬ 
presses  errors  in  force  measurement  (e.g.,  derived  from  data  such  as  those  in  Figure  3). 

Three  simulation  studies  were  performed  on  a  MicroVAX  minicomputer  using 
specially  written  software.  The  simulation  used  (pseudo-)  random  number  generators  to 
obtain  values  of  the  noise  terms.  A  sequence  of  "actual"  depths  of  cut  and  a  sequence  of 
"actual"  soil  strengths  along  a  presumed  plow  path  were  first  generated  by  Equations  43 
and  44,  respectively.  Then,  those  sequences  were  used  to  calculate  a  sequence  of  blade 
forces  using  Equation  45.  The  parameters  of  Equations  43,  44,  and  45  were  kept  the 
same  in  each  simulation.  The  sequence  of  blade  forces  was  "observed"  with  measurement 
error  (noise)  and  the  "observed"  blade  forces  put  through  a  Kalman  filter  to  simulta¬ 
neously  estimate  depth  of  cut  and  soil  strength  at  discrete  points  along  the  plow  path. 

Case  1:  Perfect  Observations 

In  Case  1  the  observations  of  depth  of  cut  (force)  were  assumed  to  be  made  without 
errors.  The  parameters  of  the  noise  term  in  Equation  45  are 

E[ Vp( t ) 1  =  0,  E(Vp(t) ]  =  0  [Eq  46] 

Figure  13  shows  the  time  stream  of  blade  force  for  the  simulated  shear  strength 
and  blade  depth  profiles  of  Figures  14  and  15,  respectively.  Figures  14  and  15  also  show 
the  estimated  shear  strength  and  blade  depths  obtained  by  applying  the  Kalman  filter  of 
Equation  27  to  the  time  stream  of  blade  force  in  Figure  13.  Scatter  plots  of  estimated 
and  true  state  variables  are  shown  in  Figures  16  and  17. 

The  accuracy  of  the  estimated  depth  (Figure  17)  was  greater  than  expected.  To 
investigate  this  finding,  force  was  plotted  against  depth  as  shown  in  Figure  18.  The 
relationship,  except  for  large  deviations  at  around  1.8  to  2.4  ft,  is  approximately  linear, 
which  explains  the  accuracy  of  the  estimated  depth.  The  reason  for  this  linear  behavior 
is  the  low  coefficient  of  variation  (10  percent)  assumed  for  undrained  shear  strength 
("Soil  Model,"  this  chapter).  Although  c  is  variable,  its  variation  is  small  compared  to  the 
mean  shear  strength  (400  psf). 

The  limiting  standard  deviation  of  depth-of-cut  estimate  (V  P  )  in  this  simulation 

co 

is  0.075  ft  (0.9  in).  This  suggests  that  depth  of  cut  might  be  controlled  to  within  an  inch 
(rtns).  This  standard  deviation,  however,  should  be  used  with  caution  since  this  value  is 
calculated  for  the  linearized  system. 
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Case  2:  Noisy  Observations  of  Blade  Depth 


Case  2  was  designed  to  assess  the  sensitivity  of  estimates  to  noise  in  the  force 
measurement.  Two  simulations  were  performed,  with  the  noise  in  the  force  model 
characterized  as 


E[vp(t)l  =  0, 


E[vp(t)]  =  40000, 


[Eq  47] 


E[vp(t)]  =  0, 


E[ v^( t ) ]  =  160000. 

r 


[Eq  48] 


Case  2a.  In  this  simulation  the  variance  of  the  white  noise  in  Equation  45  corre¬ 
sponds,  approximately,  to  a  10  percent  coefficient  of  variation  in  the  force  at  mean 
values  of  the  state  variables  (d  =  2  ft  and  c  =  400  psf). 


Figure  19  shows  the  time  stream  of  blade  force  for  the  simulated  shear  strength 
and  depth  processes  in  Figures  20  and  21,  which  are  the  same  as  for  case  1. 


Figures  22  and  23  compare  the  Kalman  filter  estimates  of  state  variables  to  true 
values.  As  expected,  the  estimates  in  case  2  are  less  accurate  than  case  1,  since  errors 
in  the  force  measurement  are  translated  into  estimation  errors  of  the  state  variables. 


The  limiting  standard  deviation  of  depth-of-cut  estimate  in  this  simulation  is  0.16 
ft,  about  twice  as  large  as  in  case  1. 


Case  2b.  The  coefficient  of  variation  in  the  force  model  in  this  simulation 
(Equation  48)  is  approximately  20  percent  at  mean  values  of  the  state  parameters.  This 
value  is  the  same  as  that  calculated  from  experimental  results  in  Figure  3. 


Figure  24  shows  the  time  stream  of  blade  force  for  the  simulated  shear  strength 
and  depth  processes  shown  in  Figures  25  and  26.  Note  that  the  true  states  in  this 
simulation  [c(t)  and  d(t)]  are  identical  to  those  in  case  2a  and  for  case  1. 


Figures  27  and  28  compare  the  Kalman  filter  estimates  of  state  variables  to  true 
values.  The  estimates  are  less  accurate  than  simulation  2a,  with  the  shear  strength 
estimate  particularly  biased.  The  bias  is  a  combined  result  of  the  linearization  and  the 
large  noise  variance. 


Case  3:  Redundant  Observations 


The  characteristics  of  the  force  measurement  in  case  3  were  assumed  identical  to 
those  in  simulation  2b;  however,  it  was  assumed  that  a  second,  independent,  but  noisy 
estimate  of  depth  could  be  made.  The  second  measurement  was  modeled  as 


y(t )  =  d(t)  +  v  (c) 

y 


[Eq  49] 


with  noise  parameters 


E[vy(t)]  =  0, 


E[v2(t)]  =  v2 
y 


[Eq  50] 


This  provides  a  two-dimensional  vector  of  measurements. 
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To  assess  the  improvements  in  the  accuracy  of  depth  estimate  as  a  result  of  the 

redundant  observation,  six  simulations  were  performed  for  different  values  of  the  error 
2 

variance,  v  ,  in  Equation  50.  Table  3  summarizes  the  results  of  the  simulation.  In  the 

table,  (a/p  )  is  the  limiting  standard  deviation  in  case  2b,  and  (a/p~  )+  is  the 
00  t  °°  t 

limiting  standard  deviation  with  redundant  observations.  The  last  column  in  the  table, 

designated  (a/p  )  ,  is  an  estimate  of  (V?  )  obtained  as 

CD  e  CD  £ 


(a/eO  +  = 


(1/P  )' 


[Eq  51] 


The  implicit  assumption  in  Equation  51  is  that  depth  of  cut  is  an  independent  process 
(rQ  =  0).  The  estimate,  however,  is  an  upper  limit  for  the  correlated  process  that  depth 
of  cut  is  in  reality. 

Figure  29  compares  (a/p^  )*  to  (a/p^  )*  for  values  of  v  in  Equation  50.  The 
comparison  is  close,  with  a  maximum  6  percent  deviation.  The  implication  of  these 
results  is  that  redundant  measurements  of  blade  depth,  even  if  noisy,  substantially 
improve  the  precision  with  which  depth  of  cut  can  be  estimated. 
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7  CONCLUSIONS 


The  principal  conclusion  of  this  study  is  that  the  depth  of  cut  for  the  COV  blade 
can  be  estimated  with  reasonable  accuracy  from  measured  blade  forces.  This  means  that 
blade  force  measurement  may  provide  a  satisfactory  approach  to  automating  blade 
control. 

In  the  accepted  equation  for  blade  force,  if  undrained  conditions  (i.e.,  rapid 
shearing)  and  a  smooth  blade  are  assumed,  the  term  for  undrained  shear  strength  of  the 
soil  being  plowed  dominates  the  calculation  of  blade  force.  Therefore,  the  effect  of 
variations  in  soil  unit  weight  can  be  and  generally  are  neglected  in  calculating  variations 
in  blade  force.  The  spatial  variation  of  undrained  shear  strength  has  been  found 
empirically  to  be  well-modeled  by  stationary,  Gaussian  processes  with  an  exponential 
convariance  function.  Thus,  soil  strengths  and  resulting  blade  loads  are  modeled  as  a 
Gauss-Markov  process.  Optimal  control  theory  was  used  to  develop  a  control  equation 
incorporating  random  variations  in  topography  and  measurement  errors,  which  could  be 
used  to  control  blade  depth  on  the  COV.  This  control  equation  requires  estimates  of 
depth  of  cut,  obtained  using  a  Kalman  filter.  Simulation  studies  indicate  that  as 
covariance  distances  in  soil  properties  and  topographic  variations  increase,  the  precision 
of  estimates  of  blade  depth  increases.  Similarly,  as  the  variance  of  observation  errors 
increases,  the  precision  of  blade  depth  estimates  decreases. 

The  simulation  studies  were  based  on  a  linearized  form  of  the  blade  force  equation 
which  was  adopted  to  simplify  the  control  equations.  The  error  introduced  by  the 
linearization  was  found  to  be  small.  Parameter  values  for  the  simulations  were  selected 
to  conform  to  operation  specifications  for  the  COV  and  to  empirically  typical  soil 
properties.  When  measurements  were  assumed  to  be  error  free,  the  standard  error  of  the 
estimate  of  depth  of  cut  was  found  to  be  as  low  as  0.9  in.,  much  smaller  than  expected. 
When  measurement  noise  was  increased  to  a  10  percent  coefficient  of  variation  (i.e.,  10 
percent  of  the  force),  the  standard  error  in  the  estimate  of  depth  of  cut  rose  to  1.9  in., 
still  considered  modest.  Further  increases  in  measurement  noise  led  to  further  increases 
in  the  standard  error  of  the  depth  of  cut  estimate. 

Simulation  experiments  were  also  undertaken  presuming  that  a  second,  independent 
measurement  of  blade  cut  could  be  made  in  conjunction  with  blade  force.  Significant 
increases  in  estimate  precision  (i.e-,  decreases  in  standard  error)  can  be  realized  using 
such  a  redundant  measurement,  even  if  the  noise  in  that  second  measurement  is  large. 

While  direct  force  measurement  appears  a  fruitful  approach  to  controlling  blade 
depth,  increased  precision  can  be  gained  by  using  one  or  more  redundant  sensor  systems. 
These  redundant  systems  need  not  be  highly  accurate  to  be  effective. 


METRIC  CONVERSIONS 

1  in.  =  2.54  cm 

1  ft  =  0.3048  m 

1  psf  =  4.882  x  10'4  kg/mZ 
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Table  1 


Contribution  to  Total  Force 


c/y  P 


0.5yd2  2cd  qd 

Lift  Shear  Surcharge 


2 

10096 

16.596 

6796 

16.596 

3 

100 

12.5 

75 

12.5 

4 

100 

10.0 

80 

10.0 

5 

100 

8.3 

83 

8.3 

6 

100 

7.0 
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Figure  3.  Measured  and  theoretical  values  of  soil  resistance 
for  bulldozer  blade,  c  =  3.8  Ib/in./in.,  y=  0.0. 
(Source:  J.  R.  O'Callaghan  and  P.  J.  McCullen, 
"Soil  Mechanics  in  Relation  to  Earth-Moving 
Machinery,"  Institution  of  Mechanical  Engineers— 
Proceedings  (Symposium  on  Earth-Moving  Machi¬ 
nery),  Vol  179,  part  3F  (1964-65),  pp  23-29. 
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Figure  23.  Estimated  vs.  true  depth  for  ease  2a 


Table  3 


Effect  of  Redundant  Observation  on  Estimation  Accuracy 

( Vp  (VT  )+  (  Vp  )+  %  Error 


•3 
*1 


APPENDIX: 


LIST  OF  VARIABLES 


N  ,  N  ,  N 
4>  c*  q 


Bcc(r) 


shear  stress  on  the  failure  surface 


cohesion  intercept;  equal  to  undrained  shear  strength  under  assumptions 
of  this  report 


normal  stress  on  the  failure  surface 


soil  friction  angle 


force  on  the  blade 


unit  weight  of  soil 


depth  of  cut 


surcharge  on  the  soil 


constants  depending  on  the  soil  properties,  blade  roughness, 
and  geometry 


covariance  function  at  separation  distance 


variance  of  the  fluctuations  around  the  mean 


distance  between  any  two  points 


scale  parameter  called  the  covariance  distance 


expectation  operator 


covariance  function  of  undrained  shear  strength  of  soil  at  separation 
distance  r 


Gaussian  white  noise 


variance  of  undrained  soil  strength 


independent  Gaussian  sequence 


vector  of  independent  identically  distributed  (iid)  measurement  errors 


e"1/ro  ;  Eq.  12 


n*n  one-step  transition  matrix;  Eq  15 


n-dimensional  state  vector 


m-dimensional  control  vector 
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n-dimensional  vector  of  random  perturbation 
independent  Gaussian  white  noise  sequence 


n*m  one-step  transfer  matrix 


arbitrary  positive  definite  matrices 


performance  criterion 


feedback  matrix,  depending  on  A,  B,  f,  and  r 


identity  matrix 


matrix  transpose 


covariance  matrix  of  the  depth-of-eut  perturbations  w(t) 


initial  covariance  matrix  of  the  state  vector  x(o) 


deterministic  matrix 


vector  of  observations  of  the  states  x(t) 


deterministic  function 


covariance  matrix  of  the  measurement  errors  v(t) 


mean  estimate  of  x(t)  conditioned  on  the  measurements 
yXO),  ^r(l),  . . .  ,/(t) 


asymptotic  variance  of  the  estimated  state  vector  for  long 
lengths  of  plowing 


asymptotic  standard  error  of  the  estimated  state  vector  for  long 
lengths  of  plowing 
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